First-principles calculations of hot-electron lifetimes in metals 
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First-principles calculations of the inelastic lifetime of low-energy electrons in Al, Mg, Be, and Cu 
are reported. Quasiparticle damping rates are evaluated from the knowledge of the electron self- 
energy, which we compute within the G W approximation of many-body theory. Inelastic lifetimes 
are then obtained along various directions of the electron wave vector, with full inclusion of the 
band structure of the solid. Average lifetimes are also reported, as a function of the electron energy. 
In Al and Mg, splitting of the band structure over the Fermi level yields electron lifetimes that 
are smaller than those of electrons in a free-electron gas. Larger lifetimes are found in Be, as a 
result of the characteristic dip that this material presents in the density of states near the Fermi 
level. In Cu, a major contribution from d electrons participating in the screening of electron-electron 
interactions yields electron lifetimes which are well above those of electrons in a free-electron gas 
with the electron density equal to that of valence (4s 1 ) electrons. 
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I. INTRODUCTION 

Low-energy excited electrons in metals, with energies 
larger than ~ 0.5 eV above the Fermi energy, experi- 
ence strong electron-electron (e-e) scattering processes. 
Although inelastic lifetimes of these so-called hot elec- 
trons have been investigated for many years on the|-ba=. 
sis of the free-electron-gas (FEG) model of the solid,EJiL£l 
time-resolved two-photon photoemission (TR-2PPE) ex- 
periments have shown the key role that baudj-atructure 
effects may play in the decay mechanism.liilc2l First- 
principles calculations of hot-electron lifetimes that fully 
include the band structure of the solid have been-reported 
only very recently for aluminum and copper .E3JE3 These 
calculations^ show that actual lifetimes are the result of 
a delicate balance between localization, density of states, 
screening, and Fermi-surface topology, even in the case 
of a free-electron-like metal such as aluminum. 

In this paper, we report first-principles calculations of 
inelastic lifetimes of excited electrons in a variety of real 
metals. We start with free-electron-like trivalent (Al) 
and divalent (Mg) metals, and then focus on divalent Be 
and the role that d electrons play in a noble metal like 
Cu. First, we expand the one-electron Bloch states in a 
plane-wave basis, and solve the Xohn-Sham equation of 
density-functional theory (DFT)E3 by invoking the local- 
density approximation (LDA) for exchange and correla- 
tion (XC). The electron-ion interaction is described by 
means of non-local, norm-conserving ionic pseudopoten- 
tials, and we use the one-electron Bloch states to eval- 
uate the screened Coulomb interaction in the random- 
phase approximation (RPA).cil We finally evaluate the 
lifetime of an excited Bloch state from the knowledge of 



the imaginary part of the electron self-energy, which we 
compute, within the GW approximation of many-body 
theory.Ej Our calculations indicate that scattering rates 
may strongly depend, for a given electron energy, on the 
direction of the wave vector of the initial state. Also, 
average lifetimes, as obtained by averaging over all wave 
vectors and bands with the same energy, are found to 
deviate considerably from those derived for a FEG. 

The rest of this paper is organized as follows: Explicit 
expressions for the electron decay rate in periodic crystals 
are derived in section II, within the GW approximation of 
many-body theory. Calculated inelastic lifetimes of hot 
electrons in Al, Mg, Be, and Cu are presented in section 
III, and the conclusions are given in section IV. Atomic 
units are used throughout, i.e., e 2 = % = m e = 1. 



II. THEORY 

Take an inhomogeneous electron system. In the frame- 
work of many-body theory,CJ the damping rate r^ 1 of an 
excited electron in the state </>i(r) with energy Ei is ob- 
tained from the knowledge of the imaginary part of the 
electron self-energy, £(r,r';i?j), as 

rr 1 = -lJdr J dr'# (r)Im £(r, r'; ^)^(r'). (1) 

In the GW approximation] 2 ! one considers only the 
first-order term in a series expansion of the self-energy in 
terms of the screened Coulomb interaction: 

£(r,r';^) = ^ / dE G{v,v' ■ E % - E)W{v,v' ■ E), (2) 
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where G(r,r';Ei — E) represents the one-particle Green 
function and W(r,r';E) is the time- ordered screened 
Coulomb interaction. After replacing the Green function 
(G) by the zero order approximation (G ), the imaginary 
part of the self-energy can be evaluated explicitly: 

ImE(r,r';£,) =^</$(r')ImW(r,r / ;w)0/(r), (3) 
/ 

where to — Ei ~ Ef represents the energy transfer, the 
sum is extended over a complete set of final states <frf(r) 
with energy Ef (Ep < Ef < Ei), Ep is the Fermi energy, 
and 

W(r, r'; u) = v(r — r') + J dri J dr 2 u(r — ri) 

xx(ri,r 2 ;w)w(r 2 - r'). (4) 

Here, v(r — r') represents the bare Coulomb interaction, 
and x(r, r';o;) is the density-density correlation function 
of the solid. 

In the framewor k <af t ime-dependent density-functional 
theory (TDDFT),blo the density-density correlation 
function satisfies the integral equation 

x( r X; w )= X°(r,r';o;) + J dr x /dr 2 x°(r,ri;u;) 

x[u(n -r 2 ) + K xc (r 1 ,r 2 ;Lu)}x(r 2 ,r';uj), (5) 

where x°( r 7 r '; w ) is the density-density correlation func- 
tion of noninteracting Kohn-Sham electrons, as described 
by the solutions of the time-dependent counterpart of 
the Kohn-Sham equation. In usual practice, these am- 
plitudes are approximated by standard LDA wave func- 
tions. The kernel K xc (ti, r 2 ; u), which accounts for the 
reduction in the e-e interaction due to the existence 
of short-range exchange-correlation (XC) effects, is ob- 
tained from the knowledge of the XC energy functional. 
In the RPA, this kernel is taken to be zero. 

For periodic crystals, one may introduce the following 
Fourier expansion for the screened interaction of Eq. fl): 



BZ 



q G,G' 
-1 



xu G (q)e GG ,(q,w), 



(6) 



where the first sum is extended over the first Bril- 
louin zone (BZ), G and G' are reciprocal lattice vec- 
tors, Q is the normalization volume, v G (q) represent the 
Fourier coefficients of the bare Coulomb interaction, and 
e G 1 G l (ij w ) are the Fourier coefficients of the inverse di- 
electric function, 

e c!G'(q^) = § G,G> +XG,G'(q,w)l>G'(q). (7) 

Within RPA, 

£G,G'(q,w) = 5 G ,G' - x GiG /(q,w)i;G'(q), (8) 



where x G g' (l) w ) are the Fourier coefficients of the 
density-density correlation function of non-interacting 
Kohn-Sham electrons (see, e.g., Ref. p^ ). 

After introduction of the Fourier representation of Eq. 
(0) into Eq. (||), and in the limit that the volume of 
the system ft becomes infinite, one finds the following 
expression for the damping rate of an electron in the 
state 4>u.,m( T ) with energy E^ ni : 



-'-is: 



G,G' 



B* f (q + G)B lf (q + G') 

|q + G| 2 



xlm -e G ^ G ,(q,w) 



and 



(9) 



where w = E^ ni - Ek-q,; 

B if (q + G) = ydr^ )ni (r)e i ^+ G )- r k _ q , n/ (r). (10) 

Couplings of the wave vector q + G to wave vectors 
q + G' with G ^ G' appear as a consequence of the 
existence of electron-density variations in real solids. If 
these terms, representing the so-called crystalline local- 
field effects, are neglected, one can write 



|B l/ (q+G)| 2 Im[e G , G (q, W ) 



|q + G| 



|e G , G (q,w)| 2 

(11) 



Within RPA, Eq. (§) yields 



Im[e G ,G(q,w)] = 



2nv G (q) 

n 



BZ 



2_j 5Z(/ k ^« _ /k+q,n') 



k n.ro' 



x|(0 k ,„|e- i ^+ G )- r |0 k+q ,„,}| 2 ( 5( W -E k+ ^ n , +E Kn ). (12) 

Hence, the imaginary part of e GG (q, u>) represents a 
measure of the number of states available for real tran- 
sitions involving a given momentum and energy trans- 
fer q + G and u>, respectively, which is rcnormalized by 
the coupling between initial and final states. The factor 
|eG,G(q, w)p 2 in Eq. ( p"l| ) accounts for the screening in 
the interaction with the probe electron. Initial and final 
states of the probe electron enter through the coefficients 
%(q+G). 

If all one-electron Bloch states entering both the coeffi- 
cients Bif(q + G) and the dielectric function e G;G /(q, u>) 
were represented by plane waves, then Eqs. (§) and @) 
would exactly coincide with the GW scattering rate of 
excited, electrons in a F-EG, as obtained by Quinn and 
Ferrelltl and by Ritchie.cl For hot electrons with energies 
very near the Fermi level (Ei w Ep)j-this result yields, 
in the high-density limit (r s _<< I),c3 the well-known 
formula of Quinn and Ferrelljj 



263 rJ 5/2 (Ei ~ E F )- 2 eV 2 fs. 



(13) 



For a detailed discussion of the range of validity of this 
approach, see Ref. G9L 
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We note that the decay r" 1 of hot electrons in pe- 
riodic crystals depends on both the wave vector k and 
the band index m of the initial Bloch state. Neverthe- 
less, we also define t~ 1 {E) 1 as the average of r _1 (k, n) 
over all wave vectors and bands lying with the same 
energy in the irreducible wedge of the Brillouin zone 
(IBZ). Decay rates of hot electrons lying outside the IBZ 
are considered by simply using the symmetry property 
r _1 (S'k, n) = r -1 (k, n), where S represents an operator 
of the point group of the crystal. 

For the evaluation of the polarizability x G G / (Q> w ) an( i 
the coefficients -E>i/(q + G), Eq. (^p|), we use the self- 
consistent LDA eigenfunctions of the one-electron Kohn- 
Sham hamiltonian of DFT, which we first expand in a 
plane wave basis, 



G 



(14) 



The electron- ion interaction is described by means, p£ 
non-local, norm-conserving ionic pseudopotentialsJ23 : Eil 
and the XC potential is obtained in the LDA with use of 
the Perdew-ZungerEa uarametrization of the XC energy 
of Ceperley and Alder £3 

Well-converged results have been found with the intro- 
duction in Eq. (14) of kinetic-energy cutoffs of 12, 6, and 
20 Ry for Al, Mg, and Be, respectively. In the case of 
Cu, all 4s 1 and 3d 10 Bloch states have been kept as va- 
lence electrons in the pseudopotential generation, and an 
energy-cutoff as large as 75 Ry has been required, thereby 
keeping ~ 900 plane waves in the expansion of Eq. (14). 
Though all-electron schemes, such as the full-potential 
linearized augmented plane- wave (LAPW) method,Ej are 
expected to be better suited for the description of the 
response of localized d electrons, the plane- wave pseu- 
dopotential approach has already been successfully in- 
corporated-jn the description of the dynamical response 
of copper Ej 

Samplings over the BZ required for the evaluation of 
both the dielectric matrix and the hot-electron decay rate 
have been performed on Monkhorst-Pack (MP) meshes:E£l 
20 x 20 x 20 for Al, 24 x 24 x 10 for Mg, 24 x 24 x 16 for Be, 
and 16 x 16 x 16 for Cu. For hot-electron energies under 
study (E - E F ~ 0.5 - 4.0 eV), the inclusion of up to 40 
bands has been required, and the sum in Eq. ( jj"jj ) has 
been extended over 15 G vectors of the reciprocal lattice, 
the magnitude of the maximum momentum transfer q + 
G being well over the upper limit of ~ 2qp (qp is the 
Fermi momentum). For the evaluation of hot-electron 
lifetimes from Eq. @, with full inclusion of crystalline 
local-field effects, dielectric matrices as large as 40 x 40 
have been considered. 



III. RESULTS AND DISCUSSION 



A. Aluminum 

Due to the free-electron-like nature of the energy 
bands of face-centered cubic (fee) aluminum [see Fig. 
1], a simple metal with no d bands, the impact of 
the band structure on the electronic excitations had 
been presumed for mar w n y ears to be small. How- 
ever, X-ray measurementsE3Ea and careful first-principles 
calculati ons . of the dynamical density-response of this 
materiaEacj have shown that band-structure effects can- 
not be neglected. Full band-structure calculations of the 
electronic stopping power of Al for slow ions have shown, 
that the energy loss is ~ 7% larger than that of a FEG.EJ 
Our present calculations indicate that actual hot-electron 
lifetimes in Al are ~ 35% smaller than those of electrons 
in a FEG. 

Our ab initio GW-RPA calculation of the average life- 
time t(E) of hot electrons in Al, as obtained from Eq. 
(^) with full inclusion of crystalline local-field effects, is 
presented in Fig. 2 by solid circles. The GW-RPA life- 
time of hot electrons in a FEG with the electron density 
equal to that of valence electrons in Al (r s = 2.07) is 
exhibited in the same figure, by a solid line. Our calcu- 
lations indicate that the lifetime of hot electrons in Al 
is, within RPA, smaller than that of electrons in a FEG 
with r s — 2.07 by a factor of ~ 0.65. We have performed 
band-structure calculations of Eq. (^|) with and without 
[see also Eq. ([Ll|)] the inclusion of crystalline local-field 
corrections, and have found that these corrections are 
negligible for electron energies under study. This is an 
expected result, since Al crystal does not present strong 
density gradients nor special electron density directions 
(bondings). 

In order to understand the origin of band-structure 
effects on hot-electron lifetimes in Al, we first focus on 
the role that the band structure plays in the creation of 
electron-hole (e-h) pairs. Hence, we evaluate hot-electron 
lifetimes from either Eq. (^) or Eq. ( pi] ) by replacing the 
electron initial and final states in \Bif(q + G)| 2 by plane 
waves (plane- wave calculation). The result we obtain 
with full inclusion of the band structure of the crystal in 
the evaluation of Ime G G /(q,w) is represented in Fig. 2 
by open triangles. Due to splitting of the band structure 
over the Fermi level, new channels are opened for e-h pair 
production, and band-structure effects tend, therefore, 
to decrease the lifetime of very-lpw-energy electrons by 
~ 5%, as in the case of slow ionsMd 

In the case of moving ions, differences between 
actual decay rates and those obtained in a FEG 
only enter through the so-called energy-loss matrix, 
Im[— e G G /(q, lj)]. However, hot-electron decay rates are 
also sensitive to the actual initial and final states entering 
the coefficients f?j/(q+ G). Differences between our full 
(solid circles) and plane-wave (open triangles) calcula- 
tions come from this sensitivity of hot-electron initial and 
final states on the band structure of the crystal, showing 
that the splitting of the band structure over the Fermi 
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level now plays a major role in lowering the hot-electron 
lifetime. 

Scaled lifetimes, t{E) x (E — Ep) 2 , of hot electrons in 
Al, as obtained from our full band-structure calculation 
(solid circles) and from the FEG model with r s = 2.07 
(solid line), arc represented in the inset of Fig. 2. In the 
limit E — > Ef, the available phase space for real tran- 
sitions is simply E — Ep, which yields the (E — Ep)~ 2 
quadratic scaling of very- low-electron energies in a FEG, 
as predicted by Eq. ( |l3| ) (dashed linc).c3 However, as 
the energy increases momentum and energy conservation 
prevents the available phase space from being as large as 
E — Ef, and the lifetime of electrons in a FEG depar- 
tures, therefore, from the (E — Ep)" 2 scaling. For ener- 
gies under study, band-structure effects in Al are found to 
be nearly energy-independent; hence, our calculated life- 
times are found to approximately scale as in the case of 
electrons in a FEG, and they slightly departure, there- 
fore, from the (E — Ep)" 2 scaling. The agreement at 
E — Ep ~ 3 eV between actual lifetimes and those pre- 
dicted by Eq. (13) is simply due to the nearly thorough 
compensation, at these energies, between the departure 
of this formula from full free-electron-gas RPA calcula- 
tions and band-structure effects. 

Although the energy bands of Al show an overall sim- 
ilarity with the fee free-electron band structure, in the 
vicinity of the W-point there are large differences be- 
tween the two cases [see Fig. 1]. At this point, the free- 
electron parabola opens an energy gap around the Fermi 
level and splits along the WX direction into two bands 
(Z\ and Z4) with energies over the Fermi level. We have 
calculated lifetimes of hot electrons in these bands, with 
the wave vector along the WX direction. The results of 
these calculations are exhibited in Fig. 2, as a function 
of energy, by dotted {Z\) and long-dashed lines (Z4). Al- 
though hot electrons in the Z\ band have one more chan- 
nel available to decay along the WX direction, the band 
gap at the M^-point around the Fermi level results in hot 
electrons living longer on the Z\ than on the Z± band. 
When the hot-electron energy is well above the Fermi 
level (E — Ep > 4eV), both Z\ and Z^ lifetimes nearly 
coincide with the average values represented by solid cir- 
cles. We have evaluated hot-electron lifetimes along other 
directions of the wave vector, and have found that dif- 
ferences between these results and average lifetimes are 
not larger than those obtained along the WX direction. 
This is in disagreement with the calculations reported 
by Schone et al.Li In particular, the bending of the hot- 
electron lifetime along the WL direction at ~ 1 eV re- 
ported in Ref. 22 is not present in our calculations. The 
origin of this discrepancy is the crossing near the Fermi 
level between bands Q+ and Q_ along the WL direc- 
tion reported in Ref. ^2[ This crosidiig-is absent in the 
present [see Fig. 1] and previousEErEjO self-consistent 
band-structure calculations, which all show that at the 
VT-point of the Al band structure the level W' 2 is below 



B. Magnesium 

In Fig. 3 we show the band structure of hexagonal 
closed-packed (hep) magnesium. There is a close resem- 
blance for energies E < Ep between this band struc- 
ture and that of free electrons, though the free-electron 
parabola now splits near the Fermi level along certain 
symmetry directions. As a result, the energy-loss func- 
tion, Im[— 6q G /(q, w)], of this material is approximately 
well described within a free-electron model, and band- 
structure effects on hot-electron lifetimes enter mainly, as 
in the case of Al, through the sensitivity of hot-electron 
initial and final states on the band structure of the crys- 
tal. 

Our ab initio calculation of the average lifetime t(E) 
of hot electrons in Mg, as obtained from Eq. (^|) with 
full inclusion of crystalline local-field effects, is presented 
in Fig. 4 by solid circles, together with the lifetime of 
hot electrons in a FEG with the electron density equal 
to that of valence electrons in Mg (r s = 2.66). As in the 
case of Al, we have found that local-field corrections are 
negligible for electron energies under study. 

Scaled lifetimes of hot electrons in Mg, as obtained 
from our full band-structure calculations (solid circles) 
and from the FEG model with r s — 2.66 (solid line), are 
represented in the inset of Fig. 4. We note that actual 
lifetimes in this material scale with energy approximately 
as in the case of electrons in a FEG, thereby slightly devi- 
ating from the (E — Ep)" 2 scaling predicted by Eq. ( p"3| ) 
(dashed line). Because of splitting of the band structure 
over the Fermi level new decay channels are opened, not 
present in the case of a FEG, and band-structure effects 
tend, therefore, to decrease the lifetime of hot electrons 
in Mg by a factor of ~ 0.75. Since the splitting of the 
band structure of Mg is not as pronounced as that of Al, 
the departure of actual lifetimes from those of electrons 
in a FEG is found to be smaller in Mg than in Al. 



C. Beryllium 

Among the so-called simple metals, with no d bands, 
beryllium presents distinctive features in that its band 
structure [see Fig. 5] exhibits the largest departure from 
free-electron behaviour. Both Be and Mg have hep crys- 
tal structure and two conduction electrons per atom. 
Nevertheless, the electronic structure of Be is qualita- 
tively different from that of Mg, the Be band gaps at 
points r, H , and L being much larger than those in Mg. 
Also, the Be band gaps are located on both sides of the 
Fermi level, and the density of states (DOS) of this ma- 
terial falls to a sharp minimum near the Fermi level. 

Our band-structure calculation of the average lifetime 
t(E) of hot electrons in Be is shown in Fig. 6 by solid 
circles, as obtained from Eq. (||). The lifetime of hot 
electrons in a FEG with the electron density equal to 
that of valence electrons in Be (r s — 1.87) is represented 



4 



in the same figure by a solid line. In the inset, the corre- 
sponding scaled calculations are plotted, together with 
the results obtained from Eq. (1^) (dashed line). It 
can be seen that large deviations from the FEG calcu- 
lation occur for electron energies near the Fermi level 
{E-E F < 3eV), especially at E - E F ~ 1.4 and 1.8 eV 
where the presence of band gaps at points T and L plays 
a key role. We note that the deep departure from free- 
electron behaviour of the beryllium DOS near the Fermi 
level tends to increase the inelastic lifetime of all excited 
Bloch states. Furthermore, actual lifetimes strongly de- 
viate from the ~ (E — Ep)~ 2 scaling predicted within 
Fermi-liquid theory. This deviation comes from the con- 
tribution to the average lifetime due to Bloch states near 
the points T and L with energies close to the energy gap. 

At the T-point, the free-electron parabola opens a wide 
energy gap around the Fermi level and splits along the 
TK and TM directions into two bands, T2/T4 and S1/S3, 
respectively. The results of our calculated lifetimes of hot 
electrons in bands T2 and T4, with the wave vector along 
the TK direction, are plotted in Fig. 7 by short-dashed 
and long-dashed lines, respectively. For comparison, the 
average lifetime of hot electrons in real beryllium and in 
a FEG with r s — 1.87 are also plotted in this figure by 
solid circles and by a solid line, respectively. At very-low 
electron energies (E — E F < 1 eV), interband transitions 
yield lifetimes of hot electrons in the T2 band that are be- 
low those of electrons in a FEG, as in the case of Al and 
Mg. However, at higher energies the coupling with lower 
lying flat bands becomes small, and lifetimes of electrons 
in this (T 2 ) band are found to be above the FEG predic- 
tion. Lifetimes of hot electrons in the T4 band are also, 
at very-low electron energies (E — E F < 1.4 eV), below 
those of electrons in a FEG. Nevertheless, at energies of 
~ 1.4 eV the presence of the band gap at T yields very 
long lifetimes, especially at the level . 

We have calculated hot-electron lifetimes in bands Si 
and S3, with the wave vector along the TM direction, 
and have found results that are similar to those obtained 
for electrons in bands T2 and T4. Calculations of the 
lifetime of hot electrons in the A2 band, with the wave 
vector along the TA direction, are also shown in Fig. 7 
(dotted line). Though this is not a flat band, the presence 
of the gap at the T-point near the Fermi level results in 
hot electrons close to the level T^ living much longer than 
in the case of a FEG. The combined contribution from 
hot electrons in bands T4, Si, and A2 is the origin of 
the enhanced average lifetime (solid circles) at E — E F ~ 
1.4 eV, which corresponds to the energy of the level 
at the T-point. 

A wide band gap is also opened at the L-point, which 
originates the enhanced average lifetime at E — E F ~ 
1.8 eV. Hence, we have plotted in Fig. 8 the results of our 
calculated lifetimes of hot electrons in the S\ band, with 
the wave vector along the LH direction (dotted line), to- 
gether with the average lifetime of hot electrons in real 
Be (solid circles) and in a FEG with r s = 1.87 (solid 
line). As in the case of hot electrons in bands T4 and Si, 



the presence of the band gap at the L-point yields a very 
long average lifetime, but now at E — E F ~ 1.8 eV. A 
similar behaviour is obtained near the L-point for elec- 
trons in the R1R3 band along the LA direction, and both 
bands, Si and R1R3, contribute to the enhanced average 
lifetime at E — E F ~ 1.8 eV. In Fig. 8 we have also 
represented calculations of the lifetime of hot electrons 
in the S[ band along the HA direction (dashed line) . At 
low energies (E — E F < 2eV), the presence of the gap 
at the ff-point leads to hot-electron lifetimes along the 
HA direction that are longer than those of electrons in a 
FEG, but departure from free-electron behaviour at the 
-ff-point is not as pronounced as at the r or L points. At 
higher energies, the S[ band shows great similarity with 
the corresponding hep free-electron band, and lifetimes 
nearly coincide, therefore, with those obtained within the 
FEG model. 



D. Copper 

Copper, the most widely studied metal by TR-2PPE, 
is a noble metal with entirely filled 3cHike bands. In 
Fig. 9 we show the energy bands of this fee crystal. We 
see a profound difference between the band structure of 
Cu and that of free electrons. Slightly below the Fermi 
level, at E — E F ~ 2eV, we have d bands capable of 
holding 10 electrons per atom, the one remaining elec- 
tron being in a free-electron-like band below and above 
the d bands. Hence, a combined description of both de- 
localized 4s 1 and localized 3d 10 electrons is needed to 
address the actual electronic response of this material. 
The results presented below have been found by keeping 
all 4s 1 and 3d 10 Bloch states as valence electrons in the 
pseudopotential generation. 

Band-structure GW-RPA calculations of the average 
lifetime t(E) of hot electrons in Cu are exhibited in Fig. 
10 by solid circles, as obtained from Eq. (^) with full 
inclusion of crystalline local-field effects. The lifetime of 
hot electrons in a FEG with the electron density equal to 
that of 4s 1 electrons in Cu (r s = 2.67) is represented by 
a solid line. These calculations indicate that the lifetime 
of hot electrons in Cu is, within RPA, larger than that 
of electrons in a FEG with r s = 2.67, this enhancement 
varying from a factor of ~ 2.5 near the Fermi level (E — 
E F = l.OeV) toafactorof- 1.5 for E-E F = 3.5eV. In 
order to investigate the role that localized d bands play in 
the decay mechanism of hot electrons, we have also used 
an ah initio pseudopotential with the 3d shell assigned to 
the core. The result of this calculation, displayed in Fig. 
10 by a dotted line, shows that it nearly coincides with 
the FEG calculation; thus, d-band states play a key role 
in the hot-electron decay. 

We have performed band-structure calculations of Eq. 
(g) with and without [see also Eq. (|TT|)] the inclusion of 
crystalline local-field corrections, and have found that 
these corrections are negligible for E — E F > 1.5 eV, 
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while for energies very near the Fermi level neglecting 
these corrections results in an overestimation of the life- 
time of less than 5%. Therefore, differences between our 
full band-structure calculations (solid circles) and FEG 
calculations (solid line) come from the actual DOS avail- 
able for real excitations, localization, additional screen- 
ing, and Fermi-surface topology. 

First of all, we focus on the role that both DOS and 
coupling between Bloch states participating in the cre- 
ation of e-h pairs, i.e., localization [see Eq. (fL2|) l play 
in the hot-electron decay mechanism. Hence, we neglect 
crystalline local-field effects and present the result of eval- 
uating hot-electron lifetimes from Eq. ( pi] ) by replacing 
initial and final states in |£?i/(q+G)| 2 by plane waves and 
the dielectric function in |e G ,G(q,w)| by that of a FEG 
with r s = 2.67. If we further replaced Im [eG.G(q, w)] by 
that of a FEG, then we would obtain the FEG calculation 
represented by a solid line. The impact of the actual DOS 
below the Fermi level may be described by simply replac- 
ing the one-electron Bloch states in Eq. ( |T2| ) by plane 
waves but keeping the actual number of states available 
for real excitations. The result of this calculations is 
represented in Fig. 10 by a dashed line. This result 
is very close to that reported by Ogawa et alf3 though 
these authors approximated the FEG dielectric function 
in |cg,g(<1) w )I 2 within the static Thomas-Fermi model. 

It is clear from Fig. 10 that the actual DOS avail- 
able for real transitions yields lifetimes that are shorter 
than those obtained in a FEG model, especially for 
E — Ep > 2 eV due to opening of the d-band scattering 
channel dominating the DOS with energies ~ 2 eV. How- 
ever, if one takes into account, within a full description 
of the band structure of the crystal in the evaluation of 
Im [eG,G(q, w)], the actual coupling between initial and 
final states avaliable for real transitions, then one obtains 
hot-electron lifetimes which lie, at very- low electron en- 
ergies (E — Ep < 2.5 eV) just above the FEG curve [see 
open circles in Fig. 1 of Ref. This enhancement of 
the lifetime, even at energies below the opening of the 
d-band scattering channel, is due to the fact that states 
just below the Fermi level have a small but significant 
d component, thus being more localized than pure sp 
states. 

The combined effect of DOS and localization, which 
enters through the imaginary part of the dielectric ma- 
trix Im [e G ,G' (l) w )]; increases the lifetime of hot elec- 
trons with e nerg ies E — Ep < 2.5 eV [see open circles in 
Fig. 1 of Ref. |2l]] . As for the departure of hot-electron ini- 
tial and final states from free-electron behaviour, entering 
through the coefficients £?j/(q + G), we have found that 
it yields hot-electron lifetimes that are strongly direc- 
tional dependent, Fermi-surface shape effects tending to 
decrease the average inelastic lifetime of very- low-energy 
electrons (E - E F < 2.5 eV) [see Ref. |j. Further- 
more, the combined effect of DOS and localization, on 
the one hand, and Fermi-surface shape effects, on the 
other hand, nearly compensate. Consequently, large dif- 



ferences between hot-electron lifetimes in real Cu and 
in a FEG with r s — 2.67 are mainly due to a major 
contribution from d electrons participating in the screen- 
ing of electron-electron interactions, which is accounted 
through the factor |e G , G (q, W )P 2 m Eq. (|Tl|). 

The Fermi surface of Cu is greatly flattened in cer- 
tain regions, showing a pronounced neck in the direc- 
tion TL. Thus, the isotropy of hot-electron lifetimes in 
a FEG disappears in this material. While flattening of 
the Fermi surface along the TK direction is found to de- 
crease the hot-electron lifetime by a factor that varies 
from ~ 15% near the Fermi level {E — Ep = 1 eV) to 
- 5% for E-Ep = 3.5 eV [see also Ref. the lifetime 
of hot electrons with the wave vector along the necks of 
the Fermi surface, in the TL direction, is found to be 
much longer than the average lifetime. We have calcu- 
lated hot-electron lifetimes in the Ai band, with the wave 
vector along the TL direction, and have found the life- 
time of hot electrons at the L\ level with E — Ep = 4.2 eV 
to be longer than the average lifetime at this energy by 
-80%. 

A comparison between our calculated hot-electron life- 
times in Cu and those determined from most recent 
TR-2PPE experiments was presented in Ref. |l]. At 
E — Ep < 2 eV, our calculations are close to lifetimes 
recently measured by Knoesel et al in the very-low en- 
ergy range.113 At larger electron energies, good agreement 
between our band-structure calculations and experiment 
is obtained for Cu(110),li3 the only surface with no band 
gap in the kn = direction. 



IV. CONCLUSIONS 

We have presented full GW-RPA band-structure cal- 
culations of the inelastic lifetime of hot electrons in Al, 
Mg, Be, and Cu, and have demonstrated that decay rates 
of low-energy excited electrons strongly depend on the 
details of the electronic band structure. Though the de- 
pendence of hot-electron lifetimes in Al and Mg on the 
direction of the wave vector has been found not to be 
large, in the case of Be and Cu hot-electron lifetimes 
have been found to be strongly directional dependent. 
Furthermore, very long lifetimes at certain points of the 
BZ in Be yield average lifetimes in this material which 
strongly deviate from the ~ (E — Ep)~ 2 scaling predicted 
within Fermi-liquid theory. 

As far as band-structure effects on hot-electron ener- 
gies and wave functions are concerned, we have found 
that both splitting of the band structure and the pres- 
ence of band gaps over the Fermi level play an important 
role in the e-e decay mechanism. In Al and Mg, splitting 
of the band structure is found to yield electron lifetimes 
that are smaller than those of electrons in a FEG. On the 
other hand, large deviations of the band structure of Be 
along certain symmetry directions from the free-electron 
model near the Fermi level result in a strong directional 
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dependence of hot-electron lifetimes in this material. 

As for the presence of band-structure effects on the 
creation of e-h pairs, there are contributions from the 
actual DOS available for real transitions, from localiza- 
tion, i.e., the actual coupling between electron and hole 
states, and from screening. The combined effect of DOS 
and localization is found not to be large, even in the case 
of a noble metal like Cu with d bands. However, large 
differences between hot-electron lifetimes in Cu and in a 
FEG with the electron density equal to that of valence 
(4s 1 ) electrons are found to be due to a major contribu- 
tion from d electrons participating in the screening of e-e 
interactions. 

Crystalline local-field corrections in these materials 
have been found to be small for hot-electron energies un- 
der study. 
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FIG. 1. Calculated band structure of Al along certain sym- 
metry directions. 

FIG. 2. Hot-electron lifetimes in Al. Solid circles represent 
our full ab initio calculation of t(E), as obtained after aver- 
aging r(k, m) of either Eq. (^) or Eq. jll| ) over wave vectors 
and over the band structure for each k. The solid line repre- 
sents the lifetime of hot electrons in a FEG with r s — 2.07, as 
obtained within the full GW-RPA. Open triangles represent 
the result obtained from Eq. ( |i"l"| ) by replacing hot-electron 
initial and final states in |Bj/(q+ G)| 2 by plane waves, but 
with full inclusion of the band structure in the evaluation 
of Im [— eQ 1 G (q, w)l . Dotted and long-dashed lines represent 
the lifetime of hot-electrons in bands Z\ and Za , respectively, 
with the wave vector along the WX direction. The inset ex- 
hibits scaled lifetimes of hot electrons in Al. Solid circles and 
the solid line represent band-structure and FEG calculations, 
respectively, both within GW-RPA. The dashed line repre- 
sents the prediction of Eq. (|l3|). 

FIG. 3. Calculated band structure of Mg along certain 
symmetry directions. 

FIG. 4. Hot-electron lifetimes in Mg. As in Fig. 2 with 
r s = 2.66. 

FIG. 5. Calculated band structure of Be along certain sym- 
metry directions. 

FIG. 6. Hot-electron lifetimes in Be. As in Fig. 2 with 
r s = 1.87. 

FIG. 7. Hot-electron lifetimes in Be. Solid circles and 
the solid line represent the same quantities as in Fig. 6. 
Short-dashed and long-dashed lines represent the lifetime of 
hot electrons in bands T2 and T4, respectively, with the wave 
vector along the TK direction. The dotted line represents the 
lifetime of hot electrons in the A2 band along the TA direc- 
tion. The open square is located at the energy of the level FJ , 
showing the beginning and the end of the A2 and T4 bands, 
respectively. 

FIG. 8. Hot-electron lifetimes in Be. Solid circles and the 
solid line represent the same quantities as in Fig. 6. Dotted 
and dashed lines represent the lifetime of hot electrons in 
bands Si and S[, respectively, with the wave vector along 
the LH and HA directions. Open squares are located at the 
energy of the levels L\ and H\. 



FIG. 9. Calculated band structure of Cu along certain sym- 
metry directions. 

FIG. 10. Hot-electron lifetimes in Cu. As in Fig. 2 with 
r s — 2.67. The dashed line represents the result of replacing 
all one-electron Bloch states by plane waves (FEG calcula- 
tion) but keeping in Eq. (ltd) the actual number of states 
available for real transitions H The dotted line represents our 
full ab initio calculation of t(E), as obtained after averaging 
r(k, m) of Eq. (|^) over wave vectors and over the band struc- 
ture for each k, but with the 3d shell assigned to the core in 
the pseudopotential generation. 
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